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Abstract 

Based on a large number of smearing steps, we classify SU(3) gauge field configurations in 
different topological sectors. For each sector we compare the exact analytical predictions for the 
microscopic Dirac operator spectrum of quenched staggered fermions. In all sectors we find perfect 
agreement with the predictions for the sector of topological charge zero, showing explicitly that the 
smallest Dirac operator eigenvalues of staggered fermions at presently realistic lattice couplings are 
insensitive to gauge field topology. On the smeared configurations, eigenvalues clearly separate 
out from the rest on configurations of topological charge v, and move towards zero in agreement 
with the index theorem. 
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1 Introduction 

The study of the relation between staggered fermions and (lattice) gauge field topology has a long 
history, beginning with the work of Smit and Vink Q. In particular, it is well known that staggered 
fermions at realistic values of the gauge field coupling (3 do not show the proper relation between 
fermion zero modes and gauge field topology as it is dictated by the index theorem in the continuum. 
This established problem with staggered fermions has recently surfaced again, as lattice gauge theory 
studies have begun to test in detail the exact analytical predictions for the microscopic Dirac operator 
spectrum of staggered fermions |3|, |^, ^. The crucial point is that the whole microscopic Dirac 
operator spectrum, and not just that pertaining to the exact zero modes, has been predicted to be 
very strongly dependent on the gauge field topology |^]. It falls into one of a set of universality 
classes that in addition depend on the gauge group, the number of fermion flavors and their color 
representation [0, ^, p^ . 

In earlier studies [§, ^, ^, ^ of the spectrum of the smallest staggered Dirac operator eigenvalues 
all gauge field configurations, irrespective of their possibly non-trivial winding numbers, were simply 
bunched together, and the Dirac operator spectrum taken with respect to this full average. This implies 
an implicit sum over all topological sectors, for which the analytical predictions are very different from 
those of sectors with fixed topological index |11|. Nevertheless, absolutely excellent agreement was 



found when comparing with just the sector of vanishing topological charge = 0. On the surface 
this may seem to be just a simple consequence of the fact that also exact zero fermion modes are 
missing in the staggered formulation. The issue is, however, much more complicated. The lattice 
studies of refs. |^, ^, ^ involved the distribution of just around 10 of the lowest Dirac operator 
eigenvalues. Only positive eigenvalues were considered, since the staggered Dirac spectrum has an 
exactly it symmetric spectrum. Of these few lowest eigenvalues, typically one would expect that up to 
2-6 were actually the "would-be" zero modes, shifted away from the origin by the staggered fermion 
artefacts. What should be the distribution of these "wrong" small eigenvalues? In the Random Matrix 
Theory formulation of the problem [^, ^ there is no answer to this question, as there is no known 
way of imposing correctly almost-zero modes in the theory. In that formulation one either has or has 
not exact zero eigenvalues. This could incorrectly lead to the conclusion that as long as staggered 
fermions do not produce exact zero modes, the distribution of the smallest Dirac operator eigenvalues 
in that formulation will exactly equal that of the v = Q sector. The argument is false, because as f3 is 
increased the appropriate number of the smallest Dirac operator eigenvalues will slowly separate out, 
and move towards the origin. As this happens, the whole microscopic Dirac operator spectrum will 
continuously shift, and in this intermediate region there will certainly no longer be agreement with 
the analytical predictions for any sectors with 7^ 0, let alone the sum over all of them. 



These considerations immediately raise the question of whether already at the /?- values considered in 
refs. 1^, ^, ^ |5| there was appreciable contamination of "wrong" non-zero eigenmodes. To settle that 
issue, we report here on a high-statistics analysis (involving around 17,000 gauge field configurations) 
of the microscopic Dirac operator spectrum of SU(3) gauge theory with quenched staggered fermions. 
Correctly classifying the gauge field configurations according to topology is not a simple task, especially 
since by default we are excluded from using fermionic methods. We have chosen to do the classification 
according to the result of measuring the naive latticized topological charge 

^ = 3^ / c^^a;Tr[F^,,Fp^]e^j,p^ (1) 

on configurations obtained after a large number of so-called APE smearing steps [|l2| , details of which 
will be given below. Once classified, we have then measured the smallest Dirac operator eigenvalues 
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on the original un-smeared gauge field configurations. We in no way claim that this is an optimal 
way in which to separate out the different topological sectors, but it should at least have quite some 
overlap with other methods. In particular, it is known, on average, to produce results very similar to 
more conventional types of semi-classical cooling [13|. 



Studying the fate of the smallest Dirac operator eigenvalues with staggered fermions is here motivated 
by the need to understand the exact analytical predictions for the microscopic Dirac operator spectrum. 
But the problem is interesting also for other reasons. For instance, one needs to know the extent to 
which staggered fermions correctly couple to gauge field topology also if one wishes to compute, for 
instance, flavor singlet pseudoscalar masses in lattice QCD. We believe that the microscopic Dirac 
operator spectrum may be an excellent tool with which to assess how close simulations with staggered 
fermions are to continuum physics related to gauge field topology and the anomaly. 

This short paper is organized as follows. In the next section we briefly describe how we implement 
the smearing technique on gauge links, and show how this leads us to a classification of all gauge 
field configurations into distinct sectors, labelled by the value of i'. We also trace the evolution of the 
smallest Dirac operator eigenvalues as a function of the number of smearing steps, and show that these 
results are completely in accord with expectations. In section 3 we compute the microscopic Dirac 
operator spectrum in each of the different gauge field sectors, and compare results to the analytical 
predictions. To zoom in on precisely the smallest eigenvalue, which in topologically non-trivial sectors 
should behave very differently from the ordinary Dirac operator eigenvalues, we also compare the 
distribution of just this smallest eigenvalue with the exact analytical predictions. Finally, section 4 
contains our conclusions. 



2 Analysis of Gauge Field Topology: APE-Smearing 

We perform the statistical analysis of the gauge field topology and the Dirac operator eigenvalues 
using a total of 17454 SU(3) pure gauge configurations with volume V = 8'^ and lattice coupling 
P = = 5.1. We generate the configurations with an update consisting of 4 microcanonical 

overrelaxation sweeps over the volume followed by one pseudo-heat-bath update. The configurations 
we analyze are separated by 20 of these compound sweeps; this guarantees that the configurations are 
effectively uncorrelated. 

Let us give some technical details about the measurements of the topological charge. First, the 
naive topological density operator (|l]), when implemented on the lattice, is not truly topological: 
it is sensitive to lattice ultraviolet modes, and, if applied to the original lattice configurations, the 
(generally non-integer) value of is completely dominated by the UV noise. However, if the fields are 
smooth enough on the lattice scale, topology can be uniquely defined. 

We make the original lattice fields smoother by applying repeated APE-smearing: during one smearing 
sweep, SU(3) link variables U^{x) are replaced by 



Uf,{x) PsU{3) 



/C/^(x)+ ^ U,{x)Uf,{x + u)Ul{x + fi) 



(2) 



where the sum over v goes over both positive and negative directions. The operator -Psu(3) projects 
the 3x3 complex matrix W^j,{x), inside the brackets in (^), to SU(3), and / is an adjustable parameter. 
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Figure 1: The topological charge of one configuration, measured during the APE-smearing. After 
initial 'noise,' the charge measurement quickly stabilizes to an almost-integer value. 

The projection is performed by maximizing Tv(c/;(x)VF^(x)) over the SU(3) group elements U'^{x). 
The smearing sweeps (consisting of a smearing of each link on the lattice) are performed up to 400 
times. 



While the APE-smearing is not particularly sensitive to the actual choice of the parameter /, we used 
/ = 7 throughout our work. This value is motivated by the comparison to the RG blocking analysis 
by De Grand et al |13]. According to their results, APE-smearing is quite effective in resolving 



instantons from quantum fluctuations, while it preserves the long-distance properties of the gauge 
field configurations much better than the standard cooling algorithms. Naturally, almost any cooling 
or smearing method will destroy instantons with a size of order of the lattice spacing; however, the 
lattice topology is not well defined at these length scales anyway. 

We calculate the field tensor F^i, in Eq. (||) at lattice point x by symmetrizing over the 'clover' of the 
4 (|U, z/)-plane plaquettes which have one corner at point x. In Fig. |l] we show a typical example of 
how the topological charge (as measured with operator (||)) develops during smearing. As expected, 
the charge measurement is very noisy when the configurations are rough, but after ~ 100 smearing 
sweeps the measurement almost always stabilizes to almost an integer value, which is then preserved 
for at least several hundreds of smearing steps. 

In Fig. ^ we show the distribution of the measured topological charge after 200 smearing steps. The 
distribution is strongly peaked near integer values, as expected, with a small downward drift of the 
peaks from exact integer values. This is obviously caused by our naive FF operator; using an improved 
operator would presumably shift the peaks toward exact integer values. Nevertheless, this slight shift 
does not pose any difficulties in assigning a topological index to the smeared configurations. Indeed, it 
is fairly obvious from our figure that there is simply a trivial "renormalization" of the naive topological 
charge, which appears to be neatly quantized, not in units of integers, but in units of around ^ 0.8. 
If corrected with such a prefactor of ^ 1/0.8, our distribution of topological charges will fall on values 
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Figure 2: The distribution of the topological charge, measured from 17454 V = 8^, f3 = 5.1 configu- 
rations. 



very close to integers, and, as seen, the distribution has the expected infinite-volume Gaussian shape. 

We have measured the topological charge after 200 and 300 smearing sweeps, and accepted only config- 
urations where the measured charge remained in the neighborhood of the same "renormalized" integer 
value 0, ±1, ±2. In this way we have rejected configurations where the measured topological charge 
changed too much to be uniquely defined. In this process we have rejected 37% of our configurations, 
which simply have been discarded in all of the subsequent analysis. 

Eigenvalues of the staggered Dirac operator 



= l)e,o+l)o.e 



(3) 



are computed using the variational Ritz functional method [14|. Here rj^{x) = (— 1)S"'<m^'' are the 
staggered phase factors. Letting e{x) = (— 1)^"^", we have explicitly indicated howi^ connects even 
sites, i.e. those with e(x) = +1, with odd (those with e(x) = —1) ones, and vice versa. We remind 
the reader that the staggered Dirac operator is antihermitian with purely imaginary eigenvalues that 
come in pairs, ±iX. The operator — Z^^ is thus hermitian and positive semi-definite, and the sign 
function e(x) defined above plays the role of 75 in the continuum (i.e., this quantity anticommutes 
with^ : {i^, e} = 0). Moreover, since — /^^ does not mix between even and odd lattice sites, it suffices 
to compute the eigenvalues, on, say, the even sublattice. In fact, if 'i/'e is a normalized eigenvector 
of — i^^ with eigenvalue A^, then ipo = j0o,e'4'e is a normalized eigenvector of —I^^ with eigenvalue 
A^, and non-zero only on odd sites. (Note that there is no difficulty with the above definition of ipo, 
since we will never encounter exact zero modes). In practice, we make use of these properties, and 
compute only the (positive) eigenvalues of — Z^^ restricted to the even sublattice, and then take the 
(positive) square root. All eigenvalues to be shown in the following thus have an equal number of 
negative companions, of the exact same magnitude. 
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We have first investigated how the lowest eigenvalues behave during the smearing process described 
previously. As shown in Fig. ^, the index theorem is seen to be valid for smooth configurations: after 
~ 100 smearing sweeps, we have exactly 4 x i/ very small eigenvalues (corresponding to 4 continuum 
flavors for each staggered fermion flavor), whereas the other eigenvalues grow larger. This implies 
that in the continuum limit the eigenvalues of the staggered Dirac operator should indeed depend on 
the topology of the gauge configuration. As expected, on the smooth configurations the counting of 
flavors is precisely as in the continuum limit: we observe 4 flavor degrees of freedom associated with 
the staggered Dirac operator. 



3 The Microscopic Dirac Operator Spectrum 



Having described our procedure for classifying our gauge field configurations according to their smeared 
topological charge u, we now proceed to measure the distribution of the smallest Dirac operator 
eigenvalues. We stress that we use the smearing technique only as a means of defining the classification 
of the original un-smeared gauge field configurations. All measurements presented below are performed 
on the original, un-smeared, gauge field configurations, after having discarded those that failed to be 
classified according to the criterion described above. 
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The spectral density of the Dirac operator is given by 

pH(A) ^ (^5(A-A„)), (4) 

n 

in each sector of topological charge v. The associated microscopic Dirac operator spectrum is defined 
by enhancing the smallest eigenvalues according to the size of the lattice space-time volume V . Let 

S = lim lim (■i/'V') (5) 

m^O V^OQ 

denote the infinite- volume chiral condensate as defined in the conventional manner. One then blows 
up the small eigenvalues by keeping C, = XLV fixed as ^ — > oo, and introduces the microscopic spectral 
density 

Pr(C) . i.'"' (^) , (6) 

which then again is measured in each topological sector. For the case at hand, the microscopic spectral 
density has been computed from both Random Matrix Theory Q and from the effective finite-volume 
partition functions |^ 10|. The results agree, and the simple analytical result for the quenched theory 



(J„(x) is the nth order Bessel function), 

pH(C) = 7rp(0)| [UCY + J.-i(C) J.+i(C)J (7) 

is exact in the limit where V ^ oo and V -C Here p{Q) is the macroscopic spectral density of 

the Dirac operator, evaluated at the origin. By the well-known Banks-Casher relation, it is related to 
the chiral condensate through 7r/?(0) = S. 

In Fig. Q we compare these analytical predictions for the sectors of = 0, 1 and 2. We did not have 
enough statistics to perform a similar analysis on what we would classify as z/ = 3 configurations; 
however, as will be evident, there was no need to do this either. Our first comment concerns statistics. 
Because the configurations with +v and —v should give rise to the same microscopic Dirac operator 
spectrum, we actually end up with better statistics for the v = ±1 configurations combined. In total, 
we had 2683 configurations labelled as = 0, 4797 configurations labelled as = ±1, and 3493 
configurations labelled as = ±2. Indeed, we observe from Fig. ^ the curious fact that the agreement 
between the analytical curve for the microscopic spectral density of the staggered Dirac operator 
actually seems to be poorer on the z/ = than on the v = ±1 configurations. We attribute this solely 
to statistical fluctuations. It is quite obvious that even configurations classified as of = ±1 and 
V = ±2 topological charge give rise to a microscopic staggered Dirac operator spectrum here which 
is indistinguishable from that of the v = Q configurations. The agreement with the exact analytical 
formula for z/ = configurations (0) is extraordinarily good on all three classes of configurations. We 
have also indicated on the figure the predictions for = 1 and \v\ = 2; there is clearly no way our 
Monte Carlo data can be compatible with these predictions. Combining all data, we obtain perfect 
agreement with the analytical v = Q prediction, with very small statistical errors. 

To focus more closely on just the smallest eigenvalue, we have also compared its distribution with the 
analytical predictions for different topological sectors |15|. Let us denote the distribution of the lowest 
(rescaled) eigenvalue in a sector of topological charge v by P^'^^C,). From the general formula of ref. 
flSll one finds for the quenched theory, with our normalization convention: 
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Figure 4: The microscopic spectral density of the Dirac operator for topological sectors v = 0, 1 and 
2, and for all of the configurations combined. The continuous line is the theoretical prediction for the 
sector z/ = 0, and the dashed lines are the predictions for = 1 and 2. 



^p(0)^[/2(C)'-/i(C)/3(C) 



-CV4 



(8) 



Checking the distribution of just the smallest eigenvalue is obviously the most sensitive test of whether 
there is any appreciable contamination of would-be zero modes in the different topological sectors. 
We show in Fig. ^ the lowest eigenvalue distributions in the three different topological sectors, and 
compare them with the analytical predictions. Clearly no deviations are seen at all from the 1^ = 
prediction, in all sectors. 



Finally, one could well ask what would happen if we instead measured the distribution of the smallest 
staggered Dirac operator eigenvalues on the smeared configurations, after removing by hand those 
eigenvalues that obviously should be classified as zero modes. All indications are of course that 
such smeared configurations should produce the correct behavior of the eigenvalues in the different 
topological sectors. The reason why we have not performed such measurements on the smoothened 
configurations is that the ensemble average is not a very meaningful concept in that case. The results 
are clearly very sensitive to where we cut the APE-smearing procedure, and even with a fixed number 
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Figure 5: The distribution of the lowest eigenvalue of the Dirac operator for topological sectors = 0, 
1, 2, and for the sum over all configurations. Continuous (dashed) lines show the theoretical prediction 
for u = = 1,2). 

of smearings it is not at all obvious that this could provide us with a sensible ensemble average. 

4 Conclusions 

At a lattice coupling of /? = 5.1 the microscopic Dirac operator spectrum of staggered fermions 
in quenched SU(3) lattice gauge theory displays no deviations at all from the analytical prediction 
of the 1^ = topological sector. Even after roughly classifying all gauge field configurations into 
different topological sectors by means of the value of FF on sufficiently smeared configurations, we 
have found no deviations from the = predictions, in any of the different sectors. In view of earlier 
results, which simply bunched all gauge field configurations together independently of their proper 
assignments of topological charge ||2|, ^, ^, |5|, and which found excellent agreement on this total sample 
of configurations, these results are not extremely surprising. It is nevertheless very disturbing that at 
the presently probed lattice couplings even gauge field configurations that clearly should be classified 
as carrying non-trivial topological charge u are not at all seen as such by staggered fermions. The 
agreement with the analytical predictions for the microscopic Dirac operator spectrum in a gauge field 
sector of topological charge zero is thus a great success of the analytical framework |8|, ^ |l^, but it 
also indicates a clear failure of the staggered fermion formulation at these lattice couplings. What it 
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means is that Monte Carlo simulations with staggered fermions at these lattice couplings are oblivious 
to net gauge field topology, an essential ingredient of the dynamics of non-Abelian gauge theories. 
Such simulations are simply not mimicking the correct path integral of the continuum theory. It is 
comforting that more sophisticated fermion formulations that correctly build in the broken/unbroken 
chiral Ward identities in sectors of fixed gauge field topology now are feasible alternatives. There 
are already clear results which show that these new fermion formulations correctly reproduce the 
analytical predictions for the microscopic Dirac operator spectrum in sectors of non-trivial gauge field 
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One surprising aspect of the present work is that it shows that the typically 2-6 would-be zero modes 
of staggered fermions at these lattice spacings behave entirely as "conventional" (small) Dirac operator 
eigenvalues on topologically trivial configurations. In particular, although they consist of up to 50% of 
all the small eigenvalues we probe here, their distributions fall entirely on top of the conventional small 
eigenvalues. Certainly, these would-be zero modes will eventually, as the lattice spacing is decreased, 
separate out, and completely distort the microscopic Dirac operator spectrum. This deformation of 
the smallest eigenvalue spectrum may already have been seen in the 2-d Schwinger model |jl^. It 
is challenging to search for the corresponding onset of correct topological properties with staggered 
fermions in this SU(3) gauge theory at smaller lattice spacings. The microscopic Dirac operator 
spectrum is an excellent tool with which to measure this in a precise and quantitative manner. 
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